Percolating cluster of center vortices and confinement 
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We study the role of percolating clusters of center vortices in configurations of an Ising gauge theory in 3D. It 
is known that low energy features of gauge theories can be described in terms of an "efi'ective string picture" , and 
that confinement properties are associated with topologically non-trivial configurations. We focus our attention 
upon percolating clusters of center vortices, and present numerical evidence for the fact that these objects play 
a preminent role in confinement phenomenon, since their removal sweeps off confinement altogether. Moreover, 
numerical simulations show that the string fiuctuations, and in particular the Luscher term, are completely 
encoded in the percolating cluster. 



1. INTRODUCTION 

Center vortices are believed to play an impor- 
tant role in confining gauge theories: they disor- 
der the gauge configurations and hence make the 
Wilson loop decay with the area law. This idea 
is quite old ||l|,D, but its investigation by lattice 
methods began recently ||,^. 

We study a very simple, though non trivial, 
gauge theory, namely Z2 gauge theory in three 
space-time dimensions: it turns out to be an 
ideal laboratory to to test the relevance of center 
vortices to the confinement mechanism; in fact, 
in this theory, the gauge group coincides with 
its center, so the identification of center vortices 
can be performed without resorting to any gauge- 
fixing procedure. Moreover, the rather small size 
of configuration space enables us to get high pre- 
cision numerical results from Monte Carlo simu- 
lations. 

We focus our attention on the effect of the size 
of clusters of center vortices on confinement: the 
hypothesis we want to test, first proposed in Ref. 
|5| , is that confinement in this model is due to the 
existence of an infinite cluster of center vortices, 
i.e. a connected component, in the graph defined 
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in the dual lattice by all the center vortices, whose 
size scales linearly with the lattice volume. Only 
such a giant component in the graph can disorder 
the gauge configurations enough to produce the 
area-law decay of the Wilson loop. 

In order to shed light on the role of large clus- 
ters of center vortices, we first modify each config- 
uration of a Monte Carlo-generated ensemble in 
the confined phase by removing the largest clus- 
ter of center vortices, and show that the ensem- 
ble thus obtained does not confine. On the other 
hand, if we remove all of the small clusters of cen- 
ter vortices, and leave the largest one only, con- 
finement is preserved, albeit with a string tension 
significantly smaller than the original one. Re- 
markably, also the string corrections due to the 
quantum fluctuations of the confining flux tube 
are preserved by this transformation. 

2. THE MODEL 

The Z2 gauge model is defined by the action 

S{P)^-PY.a^ ,aa = l[a, (1) 

□ lea 

where the sum is extended to all plaquettes of a 
cubic lattice; the fundamental degrees of freedom 
are Z2 variables cr defined on links i, and each pla- 
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quette contributes the product of its links to the 
action. The model enjoys a Z2 gauge invariance: 
any local transformation consisting in flipping all 
of the cr's defined on the links that meet at a given 
site, leaves the action invariant. 

We can construct center vortices by assigning 
a vortex line in the dual lattice to each frus- 
trated plaquette in the direct lattice. The result- 
ing graph of center vortices in the dual lattice 
is gauge-invariant, and has an even coordination 
number. 

Such graph is in general made of many con- 
nected components: the clusters of center vor- 
tices. The value of a Wilson loop in a given config- 
uration is ±1 according to the number, modulo 2, 
of vortex lines that are linked to the loop. A plau- 
sibility argument for the crucial role played by the 
infinite cluster of center vortices is the following 
(see also Ref. Q): the only clusters contribut- 
ing to the Wilson loop W{C) are those linked to 
C; so, in the limit of large enough Wilson loops, 
clusters of finite size R can contribute to W{C) 
only when they are near C, and their number 
grows linearly with the length of the loop: there- 
fore, {W{C)) decays with a perimeter law, and 
the theory is deconfined. On the other hand, if 
there exists an infinite cluster of center vortices, 
then the number of vortex lines linked to the loop 
grows with its area, and the system is in the con- 
fining regime. It is easy to verify numerically that 
in the confining phase of Z2 gauge system there is 
no ambiguity in finding a cluster of center vortices 
whose size scales linearly with the lattice volume 
for large enough lattices; conversely, in the de- 
confined phase the density of the largest cluster 
decreases rapidly with the volume. 

3. ROLE OF THE LARGE CLUSTER OF 
CENTER VORTICES 

To study the relation between the infinite clus- 
ter and the string tension value, we chose to sim- 
ulate the model at /3 = 0.74883, which is well 
inside the scaling region, and for which the value 
of the string tension is known with high precision 
|,§: a = 0.01473(10). 

First, we verified that, in this confining regime, 
in the dual lattice configurations there exists a 



large cluster of center vortices, whose size grows 
linearly with the lattice volume. 

Since VEV's of rectangular Wilson loops of 
sizes i? X T in the confining regime of a gauge 
theory are expected to behave as: 



{WiR,T)) = e 



'aRT+p(B.+T)+k 



R 



(2) 



(where the term in square brackets is associ- 
ated with quantum fluctuations of the colour flux 
tube), an efficient method to extract the 
string tension a from Wilson loop data generated 
by Monte Carlo simulations consists in consider- 
ing loops of the same perimeter, but with different 
areas; we can define the following quantity: 



acff{L,n) = ^lo! 



where: 
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The estimator doff approaches the string tension 
for large L, when it is evaluated by data obtained 
from the Monte Carlo ensemble of configurations. 

Now, if we erase the largest cluster of center 
vortices (in the dual lattice) in each of these con- 
figurations, the string tension estimator cTcs ap- 
proaches zero for large L: this means that re- 
moving the "infinite" cluster leaves us with non- 
confining configurations. 

On the other hand, if we erase all clusters but 
the largest one, we still obtain confining configu- 
rations, but the string tension evaluated on these 
new configurations turns out to be significantly 
smaller than its value in the original Monte Carlo 
configurations. Another very interesting result 
(see for details) is that the corrections due to 
the flux tube fluctuations survive the elimination 
of the small clusters. 

These results may be compared with those ob- 
tained in other gauge theories [ pHl^ . 

4. CONCLUSIONS 

Our numerical results show that: 
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Figure 1. The string tension estimator evaluated 
after erasing the largest cluster in each configu- 
ration; different symbols refer to different values 
of n. 



Figure 2. The string tension estimator evaluated 
after erasing all clusters but the largest one in 
each configuration. The solid line represents the 
string tension value in the original model. 



• The largest center vortex in the dual lat- 
tice is responsible for confinement, but the 
string tension measured from configurations 
in which all other clusters have been re- 
moved does not reproduce the full string 
tension of the original theory: small clus- 
ters of vortices are unable by themselves to 
disorder the system enough to produce con- 
finement, but they do give a finite contribu- 
tion to the string tension of the full theory. 

• Quantum fluctuations of the flux tube sur- 
vive the elimination of the small clusters: 
after deletion of all the small clusters, the 
Wilson loop VEV shows the same shape de- 
pendence as in the full theory, which can be 
explained as originating by the fluctuations 
of a free bosonic string. 
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